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'■ ■ Over the last 25 years Monte Carlo programs were being developped in Cracow in the group 

guided by Prof. Stanislaw Jadach. Many of those programs became standard in their application 

domains. In the following let us review some aspects of such projects which were probably at 

rS^ • the foundation of their success. We will concentrate on mathematical aspects of their design and 

Mh| history of their construction. It is rather difficult to cover 25 years of the research in a single 

Q^i talk. That is why, I have organized my presentation around Monte Carlo PHOTOS but stressing 

QJ ' its relation to other activities and projects often realized together with Prof. Jadach. Many of 

,^ , omitted aspects will find their way into other perentations collected in this volume. 

I will concentrate on issues related to phasespace parametrization and spin amplitudes as 
used in our Monte Carlo programs such as MUSTRAAL, TAUOLA or KKMC and their similarities and 
<*" ' differences with respect to solution used in PHOTOS. 
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00 , 1. Introduction 

o 

L^ ' One of the essential steps in the construction of any algorithm for multi-particle final states 

• ^ . is the appropriate analysis of the phase space parametrization. In the PHOTOS Monte Carlo [l] 

kS I for multi-photon production, an exact phase space parametrization is embodied in an iterative 

j^ ■ algorithm, the details of which are best described in [2]. Control of the distributions and relative 

size of sub-samples for distinct numbers of final state particles requires a precise knowledge of 

the matrix elements including virtual corrections as well. In the KKMC Monte Carlo, the phase 

space generation is different, but the necessity to control matrix elements is also essential [3, il]. 

Iterative procedures for parts of amplitudes, which are at the foundation of exponentiation 

[H [5] and structure functions O [71 El 13 [10] were exploited for the sake of use in KKMC Monte 

Carlo. In particular the description of dominantly s-channel processes e~^e~ — > feZ^eTT where, 

i-channel W-exchange diagrams with gauge boson couplings, contribute to matrix elements 

provide an interesting example [11]. These studies were motivated by practical reasons, but also 

pointed at quite astonishing properties of tree-level spin amplitudes, namely that they can be 
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separated into gauge invariant parts in a semi-automated way, easy to apply in the Kleiss-Stirling 
methods p|fT3]. 

One could ask the question whether similar techniques can be used in QCD, whether they 
are of any practical use, and in fact to which degree they were already included in previous 
publications. These questions will be discussed elsewhere |14l I15| . We will not elaborate on 
these points requiring good understanding of factorization in QCD. Instead let us point to 
old, but important for me ref. [16], where properties of factorization for cross section, visualize 
themselves in a fully differential environment, even though only for QED and at first order of 
perturbartion expansion. For the sake of caution, let us mention the existence of limitations in 
such strategies, if applied to parton shower applications beyond NLO [17j. 

Our presentation is organized as follows. In Section[2]we will discuss different aspects of phase 
space parametrization, as used in PHOTOS Monte Carlo and how it compares to other prgrams. 
Discussion of approximations necessary to construct crude distributions is started in Section [2l 
Presentation of the form of first order cross section, matrix elements and approximations which 
were essential for construction of the first version of the program is given in Section [31 With 
all material collected, we will point in Section H] to mathematical properties of elements used 
in the project, which actually made it possible, even though their documentation was never of 
high priority until now. The summary in Section [5] closes the paper. 

2. Phase space 

It is of no surprize that phase space must play a central role in preparation of the algorithm of 
any Monte Carlo based on predictions originating from field theory. That is direct consequence 
of Quantum Mechanics, basic formula for cross section consist of phase space element, matrix 
element squared and the flux factor. Over many years we were stressing, in a multitude of talks 
and papers that the control of the eventual approximations is essential. Let me recall here one 
of such S. Jadach's plots, see Fig. [TJ At that time it was an achievement [1 8 [ IB] . It required 
enormous amount of work to prepare such an organization of the phase space that would be 
exact, cover complete multibody phase space, and capable to manage highly peaked distributions 
of complex structure due to collinear and soft singularities. 

As these programs are discussed elsewhere in the proceedings, let us follow here the phase 
space organization of another program originating from S. Jadach group, that is PHOTOS Monte 
Carlql]. It is also capable of covering multibody phase space distributions without any approxi- 
mation, but contrary to KKMC/KORALZ solutions conformal symmetry of the eikonal approxima- 
tion is not used. Thanks to that, this solution is closer to iterative solution used in QCD parton 
showers, but is still relatively simple to explain and formalize. 

Let us start with the explicit expression for the parametrization of an n + 1 body phase 
space in decay of the object of four-momentum P (P^ = M^), as used in PHOTOS Monte Carlo. 
As our aim is to define iterative relations, let us denote the four momenta of the first n decay 
products as ki {i = l,n) and the last n + 1 decay product as kn+i- In our case the n + 1-th 
particle will always be the real and massless photorcl. In the later steps of our construction the 
masslessnes of photons and properties of QED matrix elements will be used. 

In the following, notation from refs. |2mi21| will be used. We will not rely on any particular 
results of these papers. We only point to other, similar options for the exact n-body phase space 
parametrizations, which are also in use. 



^ The most detailed decription of the program |19l [T] , can be found in recent ref. [2] . 

^ However the construction does not rely on a photon to be massless. In principle it can be applied to define 

other phase space relations, for example the emission of an extra massive pion or emission of a pair of heavy 

particles. 
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Fig. 1. Phase space plot for the KKMC and KORALZ Monte Carlo programs. 



The Lorentz invariant phase space is defined as follows: 
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where extra integration variables: four components oi p (compensated with 6^[p — Xli ^«)) ^® 
introduced. If further, Mi...„ (compensated with (5(p^ — M^ ^)) is introduced, the element of 
the phase space takes the form: 
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The part of the phase space Jacobian corresponding to integration over the direction and en- 
ergy of the last particle (or equivalently invariant mass Mi...„ of the remaining system of l...n 
particles) is explicitly given. We will use later in the formulas mf = kf, and analogously Mi^^^n, 



defining invariant masses of ki .. .kn systems. Tlie integration over the angles 9 and (j) is de- 
fined in the P rest-frame. The integration over the invariant mass, Mi...„, is hmited by phase 
space boundaries. Anybody familiar with the phase space parametrization as used in FOWL |22j . 
TAUOLA [21j, or many other programs will find the above explanation quite standarqj- 

The question of choice of axes with respect to which angles are defined, and order in kine- 
matical construction, is less trivial. The choice for the particular option stems from necessity 
to presample collinear singularities. It is rather well known that the choice of the reference 
directions for the parametrization of the unit sphere is free, and can be used to advantage. We 
will use related, but somewhat different freedom of choice. Instead of variables 6 (p defining 
orientation of kn+i in P rest-frame we will use angles Oi (pi orienting ki (also in P rest-frame). 
The Jacobian for this reparametrization of unit sphere equals unity. 

Formula ([2|) can be iterated and provide a parametrization of the phase space with an 
arbitrary number of final state particles. In such a case, the question of orientation of the 
frames used to define the angles and the order of Mj...„ integrations (consequently, the choice 
of limits for Mj...„ integration), becomes particularly rich. Our choice is defined in ref. |1|. We 
will not elaborate on this point here. 

If the invariant mass Mi...„ is replaced with the photon energy defined in the P rest-frame, 
k^, then the phase space formula can be written as: 
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k^dk^d cos C 



X dLipSnip -^ ki...kn), (3) 



If we would have / photons accompanying n other particles, then the factor in square brackets is 
iterated. The statistical factor jy would complete the form of the phase space parametrization, 
similar to the exponent. The last formula, supplemented with definition of frames with respect 
to which angles are defined is used to define the full kinematic configuration of the event. From 
angles and energies (k-y.) of photons and also angles, energies and masses of other decay products, 
four-momenta of all final state particles can be constructed. 

If in formula ^ instead of dLipSn{p -^ ki...kn) one would use dLipSn{P —* ki...kn) the 
tangent space would be obtained. Then kn+i photon does not affect other particles' momenta 
at all, and thus has no boundaries on energy or direction. If this formula would be iterated 
then all such photons would be independent from one another as welo Energy and momentum 
constraints on the photon(s) are introduced with the relation between tangent and real n + 1- 
body phase space. The formula defining one step in the iteration reads as follow^ 

dLipSn+liP ^ki...kn, kn+l) = dLipS+' *""^^"* X iy„"+l, 

dLips+^ *''"^""* = dk^dcosed<j) X dLipSniP -^ h-kn), 

{ki, . . . ,kn+l} = T{k^,e,<P,{h, . . . ,kn}). (4) 

The W^^^ depends on details of T, and will be thus given later in formula (jlOp . To justify (H]), we 
have to convolute formula ([2|) for LipSn+i{P — > ki...kn, ^n+i) with itself (for LipSn{p — > ki...kn)): 

dM"^ 
LipSn+i{P -^ ki...kn,kn+i) = — ^^^^Lips2{P -^ kn+iP) X Lipsnip -^ h---kn) 

zvr 

Lipsnip ^ ki...kn) = — :^''^Lips2ip^ kip') X Lipsn-iip' ^ k2...kn) (5) 

zvr 



■^ The parametrizations of such a type, use properties of the Lorentz group in an explicit manner, in particular 
measure, representations and their products. That is why, they are useful, for event building Monte Carlo 
programs in phase space constructions based on boosts and rotations. 

* Expression Q would be slightly more complicated if instead of photons a massive particle was to be added. 

^ The {fci, . . . ,k„} can be identified with the event before the radiation of kj is introduced. 



and use it also for Lipsn{P —* ki...kn): 

LipSn{P —> ki...kn) = — T^^^-Lips2(-P —^ kip) x Lipsn-iip' -^ h-'-kn)- (6) 

zvr 

Note that our tangent space of variables dk-yd cos 9d<j) is unbounded from above and the limit 
is introduced by W^^^ which is set to zero for the configuations outside the phase sace. In 
principle, we should distinguish between variables like M2...n for invariant mass oi k2 ■ ■ -kn and 
^2...n for invariant mass of A;2 ■ • • kn, but in our choice for G„, Gn+i below, M2...n = ^2...n and 
Mi,,.„ is defined anyway for the n + 1-body phase space only. 

We direct the reader to refs.|19t [T] for an alternative presentation. Let us remark that 
formula dH) is quite general, many options, motivated by the properties of the matrix elements, 
can be introduced. Generally the transformation T may differ from the choice to choice quite 
a lot. The most straightforward choice can be based on any n and n + 1 body phase space 
parametrizations using invariant masses and angles (e.g. exactly as in TAUOLA [21] formulas 11 
to 13). 

If 

and 
then 



Gn ■■ Ml^„^,9i,(l)i,Ml^„^,92,(p2,---,0n-l,4>n-l ^ ki . . . kn (7) 

Gn+1 ■ k^,0,(l),Ml^,0i,(j)i,Ml^,e2,(j)2,---,Gn-l,4>n~-l ^ ^1 • • • ^n, ^n+1 (8) 



T = Gn+i{k^,e,cP,Gn\ki, . . . ,kn)). (9) 

The ratio of the Jacobians (factors X^'"^ like in formula ([2]), etc.) form the factor WJ^~^^, which 
in our case is rather simple. 



^2(27r)3 X^/^l,ml/M^,MlJM^) ' 
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because of choice for G as explained in the Appendix of ref.[2j. Note that k^ = 2M^ " ' There 

are additional benefits from such a choice. In all relations A;2 = Lk2, ..., kn = Lkn and p' = Lp' 
common Lorentz transformation L is used. Transformation L is defined by ki,ki,p' ,p' and P; 
internal relations between four vectors k2--.kn, {k2---kn) are not needed. 
Formula (JH) can be realized algorithmically in the following way: 

1. For any point in n-body phase space (earlier generated event), described for example with 
the explicit configuration of four vectors ki...kn, coordinate variables can be calculated, 
using formula ([7]). 

2. Photon variables can be generated according to Eq. ([3]). The weight Wn~^^ has to be also 
attributed. 

3. Variables obtained in this way from the old configuration and the one of a photon can be 
used to construct the new kinematical configuration for the n + 1-body final state. The 
phase space weight, which is zero for configurations outside phase space boundaries, can 
be calculated at this point from (J4ll0p and finally combined with the matrix element. 

Here we have chosen two sub-groups of particles. The first one consisted of particle 1 alone, 
and the second, of particles 2 to n combined together. Obviously in the case of 2-body decays, 
there is not much choice when construction of the first photon is performed. 



By iteration, we can generalize formula ([4]) to the case of / photons and we write: 
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In this formula we can easily localize the tangent space for the multiple photon configu- 
ration. In this space, each photon is independent from other particles' momenta. Note that 
it is also possible to fix upper boundary on k^. arbitrary high. Photons are independent one 
from another as well. Correlations appear later, thanks to iterated transformation T. The fac- 



tors W" 
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n+j— 1 



are calculated when constraints on each consecutive photon are introduced; the 
previously constructed ones are included in the n + i — 1 systen^j. 

Of course, for the tangent space to be useful, the choice of the definition of T must be 
restricted at least by the condition {fci, • • • A;„} ^ {Aii, • • • A;„} if all A;^- ^ OlH 

It is important to realize that one has to choose matrix elements on the tangent space to 
complete the construction used in PHOTOS. The number and energies of photons will be generated 
on the tangent space first. Regularization of (at least) soft singularity must be defined. Rejection, 
and event construction, is performed with the help of formula (jl]) for each consecutive photon. 
It diminishes photon multiplicity with respect to the one defined for the tangent space. Of 
course, as rejection implements changes in phase space density, a matrix element (with virtual 
corrections) of the physical space can be introduced as well. 

The treatment of the phase space presented here lies at the heart of the construction of 
PHOTOS kinematics, and was used since its beginning. It exhausts the case when there is only 
one charged particle in final state. For multiple charged particle final states new complication 
appear, because all collinear configurations need simultaneous attention, and not only the one 
along fci direction. A presampler with multichannel generation is needed. In our case we follow 
the same method as explained in ref. [21]. 

Let us now sum the above expression over /. If we add arbitrary factors f(k-y-,0-y-,(j)^.) and 
sum over I we obtain: 
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Some parts of rhs. taken alone, give crude distribution over tangent space (orthogonal set 
of variables ki,9i,(pi). Factors / must be integrable over this tangent space and regulators of 



® Configurations of k-y. which can not be resolved are reduced to the ones with that photon dropped out. 

"^ In fact further constraints have to be fulfilled to enable presampling for the collinear singularities. 
Note that variables fc-y^ , d^^ , cf>^,^ are used at a time of the m— th step of iteration only, and are not 
needed elsewhere in construction of the physical phase space; the same is true for invariants and angles 
M2...n, S\i 4'iy ■ • ■ 1 Sn-\, fpn-i ~^ fci . . . fcn of (|7l8p . whlch are also redefined at each step of the iteration. 



singularities must be introduced. We may simply request that 
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and that sum rule originating from perturbative approach (Kinoshita-Lee-Nauenberg theorem) 
can be used to control virtual corrections; both for tangent and later also final distributions. 

At this point we already have Monte Carlo solution of PHOTOS phase space. In reality, for 
that solution to work, real emission and virtual corrections need to be calculated and their 
factorization properties must be understood. That is why, choice of / is free only in principle, 
in practice it must be synchronized with those results for the sake of program efficiency. In 
case of final state QED bremsstrahlung it is rather simple, eventual complications due to QED 
corrections to rates are of no major consequences [23] for the program construction. Non leading 
corrections appear only. 

Note that this formula is very close to other ones, used in other progams or calculations. For 
example formal solution [TOl [9] of evolution equation reads 

D{x, /3,h) = S{1 -x) + PchP{x) + ^^PlhiP X P}{x) + j.l^lhiP xPx P}{x) + ■■■ (13) 

where P{x) = 6{1 - x){\ne + 3/4) + 6(1 - j; - e)^{l + x^)/{l - x) and {P x P}{x) = 

Jq dxi /q dx2d{x - XlX2)P{xi)P{x2). 

One can easily observe, that in the LL contributing regions, the phase space Jacobian's 
as used in PHOTOS trivialize pL] and lead directly to this solution. In 1994, this solution was 
truncated to second order. It was indeed profitable that solutions for similar problems were 
available in Cracow at that time. Let us give one example [23|. In this first, on multiphoton 
Monte Carlos, paper written in 1987 by S. Jadach formula (3.1) is basically the same as tangent 
space of multi-photon PHOTOS (and not much different from D{x,l3ch) discussed just above): 

a{K) =exp('^(ln;^-l)ln^ + f In. 

En=o ^ Ul=i k<,^<K '-^sik,) . . . S{kM (14) 

The difference appears in projection from this tangent space to the physical one. Classical 
solution as proposed by Jadach, use conformal symmetry, projection from eikonal (tangent) to 
physical space is performed in one step. In PHOTOS eikonal symmetry is not used. Iterative 
projection is used instead, it is somewhat similar to the one introduced in TAUOLA [25] for 
radiative corrections in leptonic tau decays. Analogies to solutions used in QCD parton shower 
algorithms can be found. 

Very important aspect of all these solutions is that the structure of singularities is the same 
in tangent and final physical space. 

3. Matrix elements 

It is out of question, that detailed analysis of MUSTRAAL Monte Carlo [16], which was a 
consequence of accidental error in copying source code from punch cards to tape, was essential for 
the design of PHOTOS program. At that time (1983) I was forced to study MUSTRAAL line after line. 



Not only the two missing lineqj of code were found, but I have studied the matrix element and 
crude distributions in all possible details. This unintentionally collected experience combined 
with importance of QED radiative corrections in phenomenology of leptonic Z couplings at the 
time of preparation for first measurements of r polarization at LEP was few years later a starting 
point for PHOTOS. 

Let us recall the properties of the Z — > l^l^j matrix element as studied by me at that early 
time and also the approximate matrix element, which was and still is used in PHOTOS. 

Let us write the explicit form of the real-photon matrix element (separated from the phase 
space Jacobians), for the e+e" —f Z'^/7* -^ /x+/i~(7) process and as used in the standard version 
of PHOTOS (published in [WIT]): 
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where: G+ = Z(p+,g+), e_ = Z(p_, g_), 

0^ = Z(7, ;U~) is defined in (/i+, ;U~)-pair rest frame. 



(15) 



For its calculation (with respect to the Born cross-section) it is enough to know the four momenta 
of the Z and its decay products. In the presented formulae we follow the notation from refs. |231 
I16j . This expression is to be compared with the exact one, taken from ref. [16] : 
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The resulting weight is rather simple, and reads: 
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For its calculation the numerical value of the electroweak couplings of Z to fermions, as well as 
information on the state from which the Z was produced is nonetheless necessary. This seemingly 
trivial requirement puts new stress on the event record: the details of the process of the Z 
production need to be coded in the event record, then correctly deciphered by PHOTOS to calculate 
the process-dependent weight. From our experience this requirement of PHOTOS may be difficult 
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Fig. 2. The comparison [23] of the standard PHOTOS (with muhiple photon emission) and the KKMC 
generator (with second-order matrix-element and exponentiation). In the left frame the invariant mass 
of the fi~^fJ,~ pair; SDP= 0.00918. In the right frame the invariant mass of the 77 pair; SDP=0. 00268. 
The fraction of events with two hard photons was 1.2659 ± 0.0011% for KORALZ and 1.2952 ± 0.0011% 
for PHOTOS. For the definition of shape difference parameter (SDP) see ^^. 
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Fig. 3. The comparisons [53] of the improved PHOTOS (with multiple photon emission) and the KKMC 
generator ( with second order matrix element and exponentiation). In the left frame the invariant mass 
of the Ai+/.t- pair; SDP= 0.00142. In the right frame the invariant mass of the 77; SDP=0.00293. The 
fraction of events with two hard photons was 1.2659 ± 0.0011% for KGRALZ and 1.2868 ± 0.0011% for 
PHOTOS. For the definition of shape difference parameter (SDP) see [5S|. 



to accept by other users of event records. The authors of event generators often choose their 
own conventions in encoding the details of hard process such as qq -^ ngZ/j*; Z/j* — > IJ,~^IJ.~ 
into the event record. 

The NLO solution for PHOTOS, as presented in ref. [23] . -would therefore be feasible -with 
some universal, standard event record, nonetheless difficult due to practical issues of interfacing. 
One should ask the question, -what is the price related to the approximation as implemented in 
public version of PHOTOS. The results for this standard and NLO improved PHOTOS are collected 
in figures [2] and [3 As one can see, improvement due to the use of exact first order matrix 
elements is unquestionable. On the other hand, the standard, easier to use, version seem to 
be sufficient in practically all phenomenological applications as well. For the time being the 
problem of the optimal choice remains rather academic. 
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Fig. 4. Results [2 from PHOTOS, standard version, and SANC for S" -^ tt^K^{'-^) decay are superimposed 
on the consecutive plots. Standard distributions, as defined in the text and logarithmic scales are used. 
The distributions from the two programs overlap almost completely. Samples of 10^ events were used. 
The ultraviolet scale, jJ-^jy, was chosen to leave total decay width unchanged by QED. 



In ref [2], we presented similar modifications in the PHOTOS kernel for the decay of B mesons 
into a pair of scalars. As one can see from the comparison of plots in figures HI [5] and [6] the 
implementation of the exact (but scalar-QED only) kernel brings a minuscule improvement in 
the agreement betvi^een PHOTOS and the reference exact simulation of SANC [27]. In this case 
both: SANC and PHOTOS are used to simulate single photon emission. (There exists no reference 
simulation with which the multi-photon version of PHOTOS could be compared.) 

For the NLO kernel in PHOTOS the results are indistinguishable from those of SANC, even at 
statistical level of 10^ events. In this case, the price paid for improvement seems to be zero, as 
there is no need for extra information to be pumped from the event record to the calculation of 
the PHOTOS weight. Actually, the exact kernel is even simpler than the standard one. 

This high precision as documented in figs. [5] and [6] is elusive: the dependencies on the 
production process may appear if form-factors (originating from some unspecified here models) 
which have to be fitted to the data. 

From the technical side, one can interpret this excellent agreement as a strong test of nu- 
merical performance of the program. The necessary studies of the exact parametrization of 
the phase space used in PHOTOS, which will also be important for future version of PHOTOS, are 
described in detail in the journal version of ref. [2]. 
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Fig. 5. Results [2] from PHOTOS, standard version, and SAMC for ratios of the i3° — + tt K~^{^) distributions 
are presented. Differences between PHOTOS and SANC are small, but are clearly visible now 



4. Mathematical aspects of the solution 

One can ask if there is anything substantial in common in all these solutions presented in 
Section [21 and whether systematization with the help of mathematical language is worth an 
effort. Indeed, at the time of writing the first versions of the programs, which are now in a wide 
use, such considerations were of low priority. In fact to a good reason: they were expected to 
slow progress and bring little. 

At present, when multidude of different solutions is available and technical complexity of 
details dominates over main principles of construction such effort may be well motivated and 
bring useful results. 

Let us look at fig. [7| where points, curved lines and surfaces on this heuristic plot represent 
consecutive manifolds of phase spaces for n, n+1, n+2 particles. Note that the dimensionality 
of manifolds is in principle counted by number of particles times dimension of Lorentz group 
representation, minus overall energy-momentum and orientation constraints. Curvature appears 
as an ultimate expansion parameter. The crude level distribution is also defined for phase spaces 
of n, n+l, n-|-2 particles but as energy momentum constraint affects only first n particles, 
the further ones constitute fiat Carthesian sub-space. One step of the iterative projections as 
presented in Section [2] is symbolically presented in fig. [8l 

Case of QED and exponetiation of multiple photon radiation is rather simple, we do not need 
to worry about topological structure which is the same for the final (physical) phase space of 
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Fig. 6. Results [2] from PHOTOS with the exact matrix element, and SANG for ratios of the i3° — > tt K~^{'y) 
distributions. Differences between PHOTOS and SANG are below statistical error for samples of 10® events. 



multiphoton configuration and for the tangent space (constructed from eikonal phase space and 
matrix elements). The projection from tangent space to the real one is trivial (at least from the 
point of view of topological properties). In case of QCD we may expect complications, on the 
other hand hadonization models simplify the task anyway as they enforce separation of colour 
in the specific way. On the other hand it may be unhelpful for the discussion of the systematic 
errors. 

There is another mathematical concept which is worth mentioning. Thanks to infrared 
sensitive regions of n+1 body phase space we obtain, in a natural way, a triangulation line for 
this n+1 body phase space manifold. In fact, structure of such induced triangulation needs 
to be (topologically) the same for tangent and physical space, the projections must match 
these triangulations. One can realize that the language of CW complexes (known in theory of 
homotopy groups) may be useful to systematize the description and to separate it into easier to 
digest parts. 

Finally let us point to nice relation between PHOTOS algorithm for single (and fixed order) 
bremsstrahlung on one side and for the multibremsstrahlung cases. The relation is a consequence 
of the properties of the tangent spaces. It can be seen from formal expansion of Poissonian 
distribution into sum of binomial ones. In the following formula we identify coefficients of 
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Fig. 7. Symbolic reprentation of phase space with up to two extra particles. Curved surface represent 
actual phase space and the flat one tangent space. The thin bands represent configurations where only 
one extra particle is added. Point in the center configuration of the Born level. It is implicitly assumed 
that particles of soft momenta do not provide much difference with respect to configurations when they 
are absent. That is why symbolically such configuration seem to coincide. 




Fig. 8. As in figl7]this plot symbolically reprents phase space with up to two extra particles. Curved sur- 
face represent actual phase space and the cylindiric one the tangent space, where projection of kinematical 
constraint of one of its dimensions was already executed. 



binomial and poissonian distributions: p = X-, q = 1—p. Powers oip denote distinct multiplicities. 
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These somewhat unexpected numerical constants, just ratios of natural numbers, provide 
trivial example of expansion of one set of special functions into another one. The consecutive 
lines of formula (|18p correspond to expansion at respectively l''*,2"'^,3^'^ and 4*^* orders. 

5. Summary 

In this talk we have presented some principles used in Monte Carlo construction. It was a 
perfect occasion to look into history of projects, often common with Prof. Jadach's. For that 
purpose iluminating mathematical aspects of the constructions seemed to be useful. They were 
one of the cornerstones in achieving quality and robustness of the results. In the presented 
talk we have concentrated on phase space and its possible description with the help of iterative 
Monte Carlo methods. Of course, main motivation of such a systematization is to search for pro- 
totypes of algorithms to be applied e.g. in QCD. Work on matrix elements was only marginally 
mentioned here. It is only starting, but some results could have been already presented now, 
see talk by Andre van Hameren. For more, I am afraid, we need to wait for some time, even 
though some promising results are already available [28^ [29] . The next anniversary Epiphany 
conference, ten years from now, will hopefully bring some nice summary on that development. 
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